Magnetic properties of a spin-3 Chromium condensate 
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We study the ground state properties of a spin-3 Cr condensate subject to an external magnetic 
field by numerically solving the Gross-Piteavskii equations. We show that the widely adopted single- 
mode approximation is invalid under a finite magnetic field. In particular, a phase separation like 
behavior may be induced by the magnetic field. We also point out the possible origin of the phase 
separation phenomenon. 
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I. INTRODUCTION 

Since the realization of Bose-Einstein condensate of 
chromium atoms there have been considerable ex- 
perimental and theoretical efforts in exploring physical 
properties of chromium condensates. Owing to the large 
magnetic dipole moment of chromium atoms, the dipo- 
lar effects was first identified experimentally from its ex- 
pansion dynamics . More remarkably, with the precise 
control of the short-range interaction using Feshbach res- 
onance, the d-wave collapse of a pure dipolar condensate 
has been observed 3. 

In the context of spinor condensates, chromium atom 
has an electronic spin s = 3, which provides an ideal 
platform for exploring even richer quantum phases as 
compared to those offered by the spin-1 and spin-2 
atoms I, i, i, 0, 11 i, [M [ill, E Hi- To date, the- 
orists have mapped out the detailed phase diagram of 
a spin-3 chromium condensate [3, [3, [111 ■ particu- 
lar, a more exotic biaxial nematic phase was also pre- 
dicted [13]. The possible quantum phases and defects 
of spin-3 condensates were also classified based on the 
symmetry considerations [13, [3- Other work on spin- 
3 chromium condensates includes theoretically studying 
the strongly correlated states of spin-3 bosons in optical 
lattices [1^ and the Einstein-de Haas effect in chromium 
condensates [isl [20j. 

Nevertheless, all the previous work concerning the 
ground state and the magnetic properties of spin-3 Cr 
condensates has adopted the so-called single-mode ap- 
proximation (SMA), which assumes that all spin com- 
ponents share a common density profile. However, the 
studies on spin-1 case show that, for an antiferromag- 
netic spinor condensate, SMA is invalid in the presence 
of magnetic field for antiferromagnetic spin exchange in- 
teraction [2h] . One would naturally question the validity 
of SMA for spin-3 condensate since the short-range inter- 
actions involved here are more complicated than those in 
spin-1 system. 

In the present paper, we study the ground state proper- 
ties of a spin-3 chromium condensate subject to a uniform 
axial magnetic field by numerically solving the Gross- 
Pitaevskii equations. We show that even though SMA 
is still valid in the absence of an external magnetic field. 



it fails when the magnetic field is switched on. More re- 
markably, we find that when the undetermined scattering 
length corresponding to total spin zero channel falls into 
a certain region, the magnetic field may induce a phase 
separation like behavior such that the peak densities of 
certain spin components do not occur at the center of the 
trapping potential. 

This paper is organized as follows. In Sec. [Tll we in- 
troduce our model for numerical calculation. The results 
for the ground state structure of a spin-3 condensate un- 
der an external magnetic field are presented in Sec. IIIII 
Finally, we conclude in Sec. |lVl 



II. FORMULATION 

We consider a condensate of N spin s — i chromium 
atoms subject to a uniform magnetic field B = Bz,. In 
mean- field treatment, the system is described by the con- 
densate wave functions '0m (m ~ —3, —2, . . . , 3). The to- 
tal energy functional of the system, _E[0m,V'm], '^^'^ be 
decomposed into two parts E — Ea+ Ei with Eg and Ei 
being, respectively, the single-body and interaction ener- 
gies. Adopting the summation convention over repeated 
indices, the single-body energy can be expressed as 



Eo= dvt 



' 2M 



1pm.', 



(1) 



where M is the mass of the atom, the trapping poten- 
tial Voxt(r) = ^Muj']_{x'^ + y"^ -\- rf'z'^) is assumed to be 
axially symmetric with rj being the trap aspect ratio, 
s = (s^,s^,s^) arc the spin-3 matrices, 5 = 2 is the 
Lande 5-factor of ^^Cr atoms, and /is is Bohr magneton. 

The collisional interaction between two spin-3 atoms 
takes the form (1, [3] 



2s 



Mnt(r,r')=<5(r-r')^55^s, 



(2) 



5=0 



where Vs projects onto the state with total spin S and 
gs = Airh^as /M with 05=0,2,4,6 being the scattering 
lengths for the combined symmetric channel S. For 
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FIG. 1: (Color online) Left panel: phase diagram of spin-3 Cr condensate in the ao-B parameter space. The shaded region 
indicates the region where phase separation occurs (see text for details). Right panel: the main characteristics of the quantum 
phases. 



^^Cr, it was determined experimentally that ag ~ 112 a^, 
04 = 58 as, and 02 = —Tub with ub being Bohr ra- 
dius 0] , while the value of oq is unknown, and we shall 
treat it as a free parameter in the results presented below. 
Making use of the relations @ 

s 

S1-S2 = ^^[5(^+l)-24], 

s 

s 

we may replace V2, Va, and Vq in Eq. (HI by 1, Si • S2, 
and (si • 82)^, such that the interaction energy functional 
becomes 

Ei = \j drn^ [C + a\e\^ + fSTrM^ + j{s)^] , (3) 

where the interaction parameters are C — —^54 + 77-94 + 
Y5-756, 7a = go - §52 + ^34 - ^56, P = -^92 - yy54 + 
jigSe, and 7 = -^52 + Y^ff4 + -^9&- Furthermore, 
n(r) = ipmi^m is the total density, 

e(r) = i\/7(00|3m;3m')V'™Vm' 
n 

is the singlet amplitude, and 

(s)(r) = -V'mSmm'Vm' 

n 

is the density of spin. Finally, 

= 2n^^"''''*'^^ + s^s*)„„/Vm', i,j = x,y,z 

is the nematic tensor, and to obtain it, we have utilized 
the relation 

((si.S2)2)=TrAA2-i(si).(s2). 



The nematic tensor was first introduced in the liquid 
crystal physics as the order parameter Af [i^l to describe 
the orientation order of the liquid crystal molecules. 
Since J\f is Hermitian, it can be diagonalized with all 
eigenvalues Aa=i,2,3 (ordered as Ai < A2 < A3) being real 
and the corresponding principle axes being mutually 
orthogonal. Unless all three eigenvalues are equal, the 
systems with two identical eigenvalues are usually refer 
to as uniaxial nematics, while those with three unequal 
eigenvalues are biaxial ones. More importantly, Aa can 
be determined by performing Stern-Gerlach experiments 
along Ba As it can be seen from Eq. ([3]), different 

quantum phases originate from the competition of 9, (s), 
and J\f. Following the discussion of Diener and Ho [ij] , 
we shall characterize phases in a spin-3 Cr condensate us- 
ing the condensate wave functions 'tjjm , singlet amplitude 
Q, spin (s), and nematic tensor Af. 

To simplify the numerical calculations, we shall fo- 
cus on highly oblate trap geometries {ij ^ 1) such that 
the condensate can be regarded as quasi-two-dimensional 
whose motion along the 2;-axis is frozen to the ground 
state of the axial harmonic oscillator. The condensate 
wave functions can then be decomposed into 

V^™(r) - (r,/^)i/4e-''-V2^^^^(p) (4) 

with p = {x,y) and </),„ being normalized to the total 
number of atoms N, i.e., / dr(j)^^iprn — N. After inte- 
grating out the z variable, Eq gives an extra constant, 
while the interaction parameters C, a, /3, and 7 are all 
rescaled by a factor (77/27r)^/^. The mean-field wave func- 
tions {ipm} are obtained by minimizing the total energy 
functional numerically using imaginary time evolution. 
We shall focus our study on the Cr line [l3| , namely only 
the scattering length uq is allowed to changed freely, since 
experimentally, it is the most relevant case. For all results 
presented in the present work, we have chosen N = 10^, 
oj± = 2tt X 100 Hz, and rj = 10. Correspon dingly, the 
dimensionless length unit a± — y/h/{Muj±) is adopted 
in throughout this paper. 

We remark that we have neglected the magnetic dipole- 
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FIG. 2: (Color online) The typical p dependence of the densities for all spin components in the phase separation region. From 
(a) to (d), the magnetic field strengths (in units of mG) are, respectively, B — 0.0211, 0.0633, 0.0844, and 0.1689. The scattering 
length ao — 5.47as is the same for all figures. The densities of those components not shown in the figures are too small to be 
seen. 



dipole interaction energy in Eq. ([3]) for simplicity, as in 
the present work, we are concentrating on investigating 
how short-range interaction and magnetic field affect the 
ground state wave function. The ignorance of dipolar 
interaction in spinor Cr condensate was also justified in 
Ref. [14i] . Moreover, we have numerically confirmed that 
for the parameter used in this paper, the dipole-dipole in- 
teraction energy is much smaller than short-ranged spin- 
dependent interaction energy when the condensate is not 
polarized by the magnetic field. 

III. RESULTS 

Figure [T] summarizes the main results of this paper. In 
the left panel of Fig. [U we present the phase diagram 
of spin-3 Cr condensate in the uq-B parameter space, 
here we have adopted the similar notations for different 
phases as in Ref. [l3|- In the right panel, we tabulate 
the major characters of each phase. We remark that the 
resolution of phase diagram is limited by step sizes of gq 
and B when we numerically scan the parameter plane, 
therefore, it is possible that more details may emerge by 
reducing the step sizes. 

A. Condensate wave functions 

The numerical results indicate that the condensate 
wave functions can always be expressed as 

<f>m{p) = V"m(p)e'''", (5) 

where the density of mth component Umip) is an axially 
symmetric function and the corresponding phase 'dm is a 
constant independent of the spatial coordinates. In case 
the external magnetic field is completely switched off, we 
find that all wave functions (j^m (r) have the same density 
profile, indicating that the SMA is valid for spin-3 con- 
densates in the absence of magnetic field. We note that 
this conclusion also holds true for spin-1 and -2 conden- 
sates. 
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FIG. 3: (color online) The instable region (enclosed by solid 
line) of a homogeneous Cr condensate and the phase separa- 
tion region (enclosed by dashed line) of a trapped Cr conden- 
sate (Same as that in Fig. [!}. 



Once the external magnetic field is applied, SMA 
quickly becomes invalid. More remarkably, as shown 
in Fig. [21 when control parameters ag and B fall into 
the shaded region in the left panel of Fig. [l] the peak 
densities of at least one of the spin components among 
TO = —3, ±2, and do not occur at the center of the trap, 
in analogy to the phase separation in a two-component 
condensate [l^, [23|. In the absence of magnetic field, the 
system is symmetric under S0(3) rotation of the spin, 
and TO — —3 component can be populated. However, 
immediately after we switch on the magnetic field, this 
SO (3) symmetry is broken such that to = ±1 spin com- 
ponents arc highly populated under a very weak magnetic 
field. As one continuously increases the magnetic field, 
the occupation number in to = —3 component increases 
with density in the margin of the trap growing faster than 
that in the center, which induces the phase separation 
like behavior. When the population in to = —3 compo- 
nent dominates, the peak densities of all spin component 
occur at the center of the trap. We remark that similar 
behavior of the wave functions also appears outside the 
Cr line [3. 

To gain more insight into the origin of the phase sep- 
aration like behavior, we consider a homogeneous Cr 
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condensate where each spin component has already con- 
densed into the zero momentum mode. The wave func- 
tion Vm(r) for phase unseparated state is then replaced 
by a uniform c-number 

tpm = Vn^rn, (6) 

where n is a real constant and £,m are complex constants. 
The ground state can be obtained by minimizing the total 
energy E subject to the normalization condition ^*„^m — 
1. In such a way, we have reproduce the phase diagram 
in Ref. To confirm that those phases are indeed 

the ground states, we introduce a new set of variables, 
(p and Cp+ij corresponding to, respectively, the real and 
imaginary parts of the wave function as 

Cp=2(3-m) + l = R.e[^„i] and Cp=2(3-m)+2 — Im[^m]. 

We then construct the Hessian matrix H = ^ . 
For a solution to be stable, the Hessian matrix must 
be positive definite [1^. In Fig. [31 we present the un- 
stable region of a homogeneous Cr condensate on uq-B 
plane. To obtain it, we have chosen the density to be 
n = 3.3 X lO^'^cm"'^ which is the peak density of the 
trapped system in our numerical calculations. One im- 
mediately sees that the unstable region of a homogeneous 
condensate roughly agrees with the phase separation re- 
gion of the trapped system, which suggests that the pos- 
sible origin of the phase separation behavior is the insta- 
bility of the phase unseparated solution. 

We emphasize that, unlike in a binary Bose-Einstein 
condensate where the emergence of phase separation is 
determined by the strengths of intra- and inter-species 
interactions, here for a given scattering length ao, the 
phase separation like behavior is induced by the magnetic 
field. 
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FIG. 4: (Color online) The typical behaviors of |0(p)| for 
ao — — 8.27aB (a), 5.47a_B (b), and 12.35as (c). In descend- 
ing order of central value, the lines in (a) correspond to the 
magnetic field (in units of mG) B = 0, 0.0244, 0.1689, 0.2533, 
and 0.3377; those in (b) correspond to B = 0, 0.0422, 0.0844, 
0.1266, and 0.19; and finally, those in (c) correspond to B = 0, 
0.0211, 0.0422, and 0.0633. 



B. Singlet amplitude 

Since the spatial independence of Q{p) is a necessary 
condition for SMA, it can also be used as a criterion to 
check the validity of SMA. As shown in Fig. [H \<d\ is a 
constant when B = 0; while immediately after the mag- 
netic field is turned on, |0| becomes spatially dependent. 
In addition, the peak value of |8| decreases continuously 
as one increases the magnetic field until it completely 
vanish. 

In Fig. m (a), |6(p)| becomes zero only after the con- 
densate is completely polarized, while in (b) and (c), it 
vanishes once the system enters the Hi phase. There- 
fore, using singlet amplitude, we may map out the phase 
boundaries between Ai and FF, Z and Hi, and Bi and 
Hi. However, O alone is incapable of determining other 
phase boundaries. Finally, we note that, for ao > S.Oas, 
the value of |0| drops much faster with the increasing 
magnetic field than that corresponding to oq < S.Oas, 
as shown below this behavior has a direct impact on the 
magnetization curve of the system. 
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FIG. 5: (Color online) The field dependence of the magne- 
tization for ao — — 59.82a_B (solid line), — 25.45as (dashed 
line), 5.47as (dash-dotted line), and 12.35as (dotted fine). 
For ao > 8.9as {a > 0), the Ai{p) curves corresponding to 
difi'erent ao's are indistinguishable. 
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C. Magnetization 

We now turn to study magnetic field dependence of the 
total magnetization. To this end, we define the reduced 
magnetization as 



M=N- 



dp{i 



(7) 



Unlike in Ref. [T^| where the total magnetization is 
conserved, here we allow it to change freely. There- 
fore, the transverse components of the spin, (s^) and 
(s^), are always zero. Figure [5] shows the field depen- 
dence of the reduced magnetization, which approaches 
—3 when B reaches the saturation field. We note that, 
for uq < 8.9aB, the behavior oiA4{B) slightly depends on 
the value of aq; while for oq > S.Oas, the magnetization 
curves corresponding to different oq's become indistin- 
guishable. Consequently, as shown in Fig. [H the satura- 
tion field for the former case is a decreasing function of 
ao, while for the latter one, it becomes a constant. The oq 
independence of the magnetization for ao > 8.9aB case 
can be qualitatively understood as follows. The scatter- 
ing length ao only contributes to the total energy through 
singlet amplitude 6. As shown in Fig. 21 for oq > 8.9aB, 
Q vanishes quickly as one increases the magnetic field, 
such that varying ao only yields a negligible effect on 
magnetization curve. With the help magnetization, we 
can further identify the phase boundary between Hi and 
FF phases. 



D. Nematic tensor 

To determine other phase boundaries, we have to rely 
on the nematic tensor. In Fig. [51 we plot typical behavior 
of the eigenvalues of nematic tensor under different ao 
and magnetic fields. As a consequence of the failure of 
SMA, Aa's are generally spatially dependent. However, 
some of their characteristics obtained under SMA remain 
to be true. 

The full ferromagnet (FF) phase occurs when the mag- 
netic field exceeds the saturation field such that only m = 
—3 component is occupied. The nematic tensor takes a 
diagonal form with Ai = A2 = | and A3 — 9. For the 
wave functions of Ai and Hi phases, only two spin states 
are populated: other than a common m = —3 state, 
TO = 3 and 2 are also occupied for, respectively, Ai and 
Hi phases. One can easily deduce that the nematic ten- 
sors of phases Ai and Hi are, respectively, diag{|, |, 9} 



and n ^diag{^ 



i + 4^2, 9n_3 + 4n2}. Since 



for both cases, Afzz is the largest eigenvalue, we have 
63 II z. Moreover, ei and 62 are spatially independent. 

For Gi phase, to ~ ±2 and states are un- 
occupied, consequently, Af becomes a block diago- 
nal matrix with one of the eigenvalues being Afzz = 
[9(713 -|- n-3) + ni + n_i]. Furthermore, from nu- 
merical results, we find that Afzz is always the smaller 
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FIG. 6: (Color online) The spatial dependence of Ai (dash- 
dotted lines), A2 (dashed lines), and A3 (solid lines) for ao = 
5A7aB (left panels) and 12.35as (right panels). The magnetic 
field strength (in units of mG) is denoted in each figure. 



eigenvalue and A2(p) 7^ A3(p), which suggests that Gi is 
a biaxial nematic phase with ei || z. 

All spin components of Bi and Z phases are popu- 
lated, they are all are biaxial nematic with three unequal 
and spatially dependent A's. The principal axes for these 
phases are also spatially dependents. The only difference 
is that for Bi phase we can identify that either ei or 62 
is perpendicular to z-axis. 



IV. CONCLUSIONS 



To conclude, we have mapped out the phase diagram of 
a spin-3 Cr condensate subject to an external magnetic 
field based on the numerical calculation of the ground 
state wave function. In particular, we show that SMA be- 
comes invalid for Cr condensates under a finite magnetic 
field. More remarkably, if the unknown scattering length 
ao falls into the region [—0.37, 8.2]aB, a phase separation 
like behavior may be induced by the magnetic field. We 
also point out that such behavior might originate from 
the instability of a phase unseparated solution. As a 
future work, we shall investigate the ground state struc- 
ture of a spin-3 Cr condensate by including the magnetic 
dipole-dipole interaction. 
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